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Abstract. - The nature of derealization in a ID system ruled by a tight-binding Hamiltonian 
is investigated. Using a local evaluation of the ground state energy, it is shown that the range of 
the derealization effects is rather limited. The method is applied to 1/2 and 1/4 filled systems. 



Introduction. - For the last few decades, physicists and chemists have put a lot of effort 
into understanding the properties of low-dimensional compounds [1-3]. A major achievement 
was the elucidation of the physical properties of polyacetylene the simplest linear polymer [4] . 
Since then, a large number of inorganic chains have been subject to extensive studies. Usually 
the chemically relevant and extremely simple model Hamiltonian, the Hiickel Hamiltonian, 
has been called for [5-7] . The use of more sophisticated Hamiltonians such as the Hubbard 
Hamiltonian led to fundamental issues since exact solutions are known for one dimensional 
systems only [8]. The traditional solid-state picture gives a privilege to derealization, pre- 
sented as a collective phenomenon. However, it is known from the solid state and the surface 
bonding that the local interactions (i.e. bonds) play a dominant role and contain a large 
part of the chemical actions [6,9]. Following this trend, we have recently demonstrated that 
a localized description of the physics of a system ruled by very short-range interactions (i.e. 
Hiickel Hamiltonian) is relevant [10,11]. Indeed, starting from strongly localized ./V-electron 
functions built from either pure Atomic Orbitals (AOs) or fully localized bond Molecular 
Orbitals (MOs), it is possible to get very good evaluations of the ground state energy of a pe- 
riodic lattice without any computational effort. The method consists in a self-consistent (SC) 
perturbation of a zeroth-order wave-function which incorporates high order effects and offers 
particular convergence properties. The charge transfer contributions are explicitly evaluated, 
a noticeable difference with the original moments method [13] which focuses on the charge 
densities. This modified coupled-cluster approach has shed the light on the crucial importance 
of the local effects which appear to be the quantitatively dominating driving forces responsible 
for the bond alternation in 1/2 filled conjugated hydrocarbons (i.e. Peierls distortion [12]). 
The short and medium-range character of derealization in the ground state clearly appeared 
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to be dominant. However, our evaluations based either on pure AOs or bond MOs did not al- 
low one to perfectly reproduce the reference tight-binding calculations. Thus, to which extent 
derealization has to be incorporated remains a challenging issue in the study of delocalizcd 
systems. Besides, the accurate determination of the low-energy physics of such systems by 
rather local methods is tempting since it opens the possiblity to determine the ground state 
energy of systems ruled by a variety of model Hamiltonians with an efficient and unexpensive 
method. 

In this paper, we shall first briefly recall the method which has been recently proposed 
by introducing the so-called set of self-consistent coupled equations (SCPEs) [10]. We will 
then follow two different strategies based on this method. The first one consists in using MOs 
delocalized over four-site units, i.e. butadiene MOs, to build the zeroth-order wave-function. 
The second aims to extend the variational space which had been limited to single charge 
transfer amplitudes [10, 11]. Higher order terms were incoporated by means of perturbation 
theory. Thus, starting from a zeroth-order wave-function based either on AOs or MOs, single 
and double charge transfers as well as second-nearest neighbour charge transfers coefficients 
will be treated self-consistcntly. For a strictly non- interacting electron system (i.e. Huckel 
Hamiltonian) , our goal is i) to look into the benefit of both approaches in the evaluation of 
the ground state energy of a 1/2-filled homogeneous chain, ii) to extract an estimation of the 
relevant derealization range, and iii) to investigate the reliability of such local evaluations for 
1/4-filled systems. We believe that in the study of periodic lattices, local pictures are very 
valuable since the information they convey is very instructive in the undertanding of bond 
formation analysis. 

Self- Consistent Perturbation Method. Let us consider the Huckel Hamiltonian H = 
X^( a 2i-i a 2i + aL a 2«-i) ruling an infinite chain of singly occupied valence orbitals. The 
nearest-neighbour interaction is determined by the hopping integral t. The derivation of the 
Huckel energy for such 1/2-filled system is very standard and leads to 4/ it w 1.273 in the unit 
of t [7] . The basic idea behind the SCPEs method is to set the electrons in strongly localized 
monoelectronic functions (fi to build a A^-electron zeroth-order determinant <E>o- Since the a 
and (3 electrons move independently, one may pair the electrons on the odd sites to give rise 
to an Atomic Charge Density Wave (ACDW) $o = ^n-Wn-x |- One may as well start from 
the bond MOs to describe a zeroth-order Bond Density Wave (BDW) determinant. Using the 
so-called intermediate normalization, the wave-function is expanded in terms of the charge 
transfer determinants 



As in Ref. [10], we define the first generation subspace Si as the determinants {^k}k^o 
which interact with <& (i.e. ($>o\H\&k) ^ 0). Such determinants {$>k}k^o are obtained 
through charge transfers T£ between adjacent units. One can easily show that the cohesive 
energy E co h is completely determined by the knowledge of the S\ determinants amplitudes, 
E co h = H o + J2i^o c iHoi- Let us briefly recall the derivation of the SCPEs starting from the 
eigen-equation for the $^ determinant : 




(1) 



fc#0 



(H kk - E)ck + H k Q + Hkici + ^ H ka c, 



'ft: 



= 



(2) 



l^k a 



where the determinants a belong to the higher generation subspaces. If the <J/ expansion is 
limited to ^-determinants, the {c a } coefficients are determined by means of second-order 
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perturbation theory. The $ Q can be generated from by either the action of a similar T} 
operation or an operation w m which does not act on <I>o. The amplitude of the former may 
depend not only on the c k ci product associated with the T ; ' TjJ parent action on $o but also 
on the c m c„ products such that T/T^ = T\T^ n . Let us mention that a major difference with 
the traditional coupled-cluster approach [14] lies in the participation of the operators, 
the role of wich has been recently investigated by means of a block-correlated coupled-cluster 
formalism [15]. The derivation of the second-order amplitudes has been given previously [10] 
and leads to : 

AE' k + AE[ x ^ AE' m + AE' n 

c a -c k c t + ^ c ™ c « AE',,, {) 

k + l (m,n)^(k,l) k + l 

t^tI=tIt^ 

where the excitation energies AE k = H kk — Hoo are shifted by the EPV {Exclusion Principle 
Violating [16]) correction to deal with higher order corrections [17] : 

AE' k = AE k - EPV(k) = AE k - ^ H ol a (4) 

i,r/* fc =o 

Let us mention that for $ a determinants being generated by a single pair T}t1 the additivity in 
the excitation energies reduces c a to the simple product c k ci, resulting in a perfect cancelation 
of the c k ciHoi term in eq. |21 Such situation is likely to occur for remote excitations T£ and 
7j , as expected from simple physical considerations. 

The i?J n T j J$ determinants amplitudes can be immediately evaluated from second-order 
perturbation theory : 



c f = E-^ c ™ ( 5 ) 



Finally, the SCPEs in terms of the {c k } k ^o are derived from eq.|21by expliciting the coefficients 
c a using eqs.|21and|5] Let us mention that the dimension of the set of coupled equations is 
equal to the number of charge transfer types in the system. 

Part of the present work consists in a SC evaluation of the previous c a and cp coefficients. 

Extension of the Method. - For a homogeneous periodic lattice, the cohesive energy E co h 
has been previously evaluated with absolute relative errors 2.6% and 1.5% using AOs or bond 
MOs, respectively [10]. Thus, it was rather challenging to start with larger units built on 
the butadiene MOs (four-site MOs) to investigate the importance of the collective effects. In 
order to clarify the presentation, the strategy will be referred to as l,L where I stands for 
the size of the building blocks (1, 2 or 4) and L for the extent of the SC treatment. L = II 
indicates that the Si and S2 non-factorizable determinants amplitudes are treated in SC way. 

In using four-site MOs, the amount of delocalization in the zeroth-order wave-function 
$0 is enhanced as compared to bond MOs. Thus, one may expect to improve the cohesive 
energy evaluation. The derivation of the Hiickel energies {e n } n =i-4 and MOs coefficients 
of butadiene is very standard [7]. Let us write {i n }n=i-4 the MOs localized on the i th 
tetrameric unit associated with energies {e n }„ = i_4. The energy per atom corresponding to 
the zeroth-order determinant <&o = |*i*i*2*2 1 reads Eo/t = (2ei + 2e2)/4i = y/E/2 w 1.118. 
Starting from $0, one can generate the Si-subspace. Since the left and right charge transfer 
absolute amplitudes between adjacent tetrameric units are equal, only four types are to be 
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considered, namely c% = c^—^, c 2 = Ci 1 — > j i , Cj = Cj 2 _>j 3 and c' 2 — Ci 2 ^j i . As mentionned 
before, the a and f3 electrons can be treated independently in the Hiickel framework (i.e. 
E a = Ep = \/2E co h). Therefore, the cohesive energy reads : 

E coh =E Q + (A Cl h + Ac 2 t 2 + Ac[t[ + 4c' 2 4)/4 

Using the expressions of the butadiene MOs in terms of the AOs, one immediately finds the 
coupling matrix elements t\ — —t' 2 — 0.224i, t[ = — 0.362i and t 2 — 0. 138i. The generation of 
the SVsubspace is straightforward but somehow lengthy and will not be fully detailed here. 
However, the three types of determinants belonging to 62 are shown in fig. ^ The first one 




i j k I i j k h i j k 

(a) (b) (c) 

Fig. 1 - Second-order determinants : (a) T^T}^, (b) with two parents, T^T^ = T^T^, 

(c) RlnTj ■ the dotted arrows indicate the different possible propagations. The black circles stand 
for the a electrons. 



corresponds to a double charge transfer generated by a single pair {T fc T ,T/}. Conversely, the 
second exemplifies the multiple generation of a double charge transfer (see eq.|3J). Finally, the 
third one gives rise to contributions evaluated by eq. The derivation of the SCPEs for a 
1 / 4- filled system follows exactly the same pattern by taking the inventory of the different S\ 
and SVsubspaces determinants. For both fillings, the determination of the nearest-neighbour 
charge transfer amplitudes implies up to four butadiene units, that is 26 AOs, without elimi- 
nating the translational property of the chain. 

Instead of increasing the degrees of freedom one can extend the SC treatment. Thus, the 
second strategy we considered is based on an enlargement of the variational space by treating 
the S± and non-factorizable S2 determinants amplitudes c and d in a SC way. Third-order 
terms are then evaluated from perturbation theory. As a zeroth-order wave-function, we first 
used the strongly localized determinant built on the AOs. The ^-subspace consists of the 
double charge transfers which cannot be factorized as product of simple charge transfers 
(see fig. EJ). 




Fig. 2 - Charge transfers generation using the AO based zeroth-order determinant. Dotted lines 
indicate the electron propagations. 



The eigen-equation for the Si determinants $i is easily derived 

(-3ct)c + t + 2t(d ~ c 2 ) = 



(0) 
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Starting from the S2 determinants the adjacent excitation gives rise to a charge 

transfer <E>^ (see fig- EJ with amplitude eg. A third-order expansion gives : 

c p = -d [H kf 3/AE' 0f3 \ = cd [3ct/7ct] 

However, can be generated by two other routes. The first one which has the same amplitude 
cp results from the action of T/ on <F,+fc- One can also start from the strictly additive double 
excitation T^T"/ (i.e. c 2 amplitude) and then propagate an electron along the j bond through 
Tj . The resulting contribution is (3ct/7ct)c 3 . Then, one has to consider the action of on 
(eq. 0) which generates a third-order determinant $ 7 (see fig. 0) with amplitude dt/4ct 
(eq. 03). Remembering that is directly coupled to the S\ determinants $i and the 

second eigen-equation reads : 

(-5ct)d + 2ct + 2t [(3ct/7ct)(2cd + c 3 ) - cd)] + 

t [dt/Acb] = (7) 

whilst the cohesive energy 2ct is fully determined by c. In the following, eqs. and [7| are 
referred to as extended SCPEs. 

In order to get a clearer understanding of the role of the operators, one can finally start 
from a bond MOs zeroth-order wave-function $0 and use a similar second generation logic 
to evaluate the cohesive energy (1 + c)t. It has been previously observed that the deviation 
from non additivity is negligible. Thus, the double charge transfer ^1* amplitudes (see 
fig. 01 can be approximated by c 2 . However, the W m operators give rise to second- nearest 
neighbour charge transfer $ik* (see fig. OJ) the amplitudes of which d are self-consistently 
calculated. Along this framework, the third-nearest neighbour charge transfer amplitudes 
are estimated from perturbation theory. Such strategy combines the use of intermediate 
fragments size (bond MOs) and an extension of the SC treatment. The extended SCPEs are 
very easily derived in terms of the first and second nearest-neighbour charge transfers c and 
d, respectively. 



Fig. 3 - Charge transfers generation using the bond MOs based zeroth-order determinant. 

For both strategies, the non- linear set of coupled equations is solved without any compu- 
tational effort. The cohesive energies for 1/2-filled systems are summarized in Table H] These 

(n) 

energies are compared to the ones of finite size clusters E co ^ the extensions of which n are 
measured by the number of AOs implied in our approach. 

First, using the four-site based determinant, our local evaluation is in excellent agreement 
with the reference value since it deviates by less than 0.02%. The extended SCPEs strategy 
starting from the AO determinant (i.e. (l,L) = (1,11)) leads to a poorer , still satisfactory, 
estimation of the cohesive energy. One should note that the delocalization effects are smaller 
since the local evaluation is limited to 6 lattice parameters. Conversely, the combination of 
fragment of intermediate size (bond MOs) and an extended SC treatment (i.e. (I, L) — (2, II)) 
displays excellent agreement with the exact value and rather inexpensive cost. Therefore, 
considering the perturbative treatment which incorporates high-order effects through the EPV 
corrections, we may conclude that the delocalization effects are effective up to roughly 10 
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Number of sites 


1 


2 4 


Extent of the / 


1.240 a ( 1.205 (10) ) 


1.253" (1.224 (14) ) 1.273 (1.246 (26) ) 


SC-treatment 


1.267 ( 1.216 (12) ) 


1.272 (1.234 (18) ) 



Table I - Cohesive energy evaluations for a 1/2-filled system (E exact — 1.273). The exponents 
indicate the number of AOs used in the finite size clusters calculations given in parenthesis. a Ref. [10]. 



lattice parameters. The benefit of our approach which explicitely includes the translational 
symmetry is significant as seen from the finite size clusters calculations (Table QJ. Second, 
inspection of the different contributions to the cohesive energy shows that a major part can 
be attributed to the excitations which are not active in a strict coupled-cluster formalism. 
Indeed, if one sets these terms to zero in the SCPEs, the relative error deteriorates (0.02% — > 
1.6%). The comparison between the SCPEs (2,1) and the extended SCPEs (1,11) strategies 
is very instructive. Indeed, the numbers of AOs implied are very similar. However, the 
use of a larger variational space in the extended SCPEs noticeably enhances the accuracy 
(1.5% — > 0.5%). Interestingly, the d/c 2 ratio is w 1.253 whereas based on the bond MOs, 
the deviation from additivity of the local charge transfers is measured by the EPV corrected 
excitation energies ratio, that is (3ct+3ct)/5ct — 1.2 [10]. Thus, using the extended SCPEs the 
interference between local excitations is somehow increased and a better estimation results. 
Similarly, inspection of the second-nearest neighbor charge transfer amplitudes shows that 
their contribution is significantly enhanced (d = 0.155) as compared to the original first 
generation approach (d = 2c(t/2)/(2t + 3ct/2) = c/(2 + 3 /2c) = 0.106) [10]. 

A very similar approach was applied to a 1/4-filled system. Starting from two-site or 
four-site MOs, the local evaluations deviate from the exact value (0.900) by 2.9% (0.874) 
and 0.55% (0.895), respectively. Even though a relative agreement with the reference exact 
value is reached , the errors are somehow larger. From band structure calculations, it is 
known that the higher the energy in the dispersion curve, the more numerous the nodes in 
the orbital. Thus, we may argue that 1/4-filled systems are more delocalized than 1/2-filled 
ones. The accuracy of our approach which introduces only short-range derealization effects 
clearly follows the derealization nature of the system. However, the extended SCPEs method 
still leads to a more accurate estimation of the 1/4-filled system energy (0.895) than the one 
obtained from a 26 AO cluster (0.831). 

Conclusions. - In summary, the method we used allowed us to investigate the dereal- 
ization effects in a ID-chain. We showed that the latter are limited since excellent evaluations 
for both 1/4 and 1/2 filled systems can be obtained from rather local evaluations. Besides, 
the method offers a way to estimate the derealization range by step by step incorporating 
the short and medium range contributions. Beyond roughly a 10-lattice parameter distance, 
their importance is vanishingly small. From a conceptual point a view, the original SCPEs 
approach [10] can be developed at will by including higher order terms in the variational 
space. Such strategy is very powerful since the benefit in terms of accuracy is significant 
whilst the SCPEs are easily derived. In the view of our results, we believe that the descrip- 
tion of important heteroatomic (A A') chains which are more localized by essence would not 
be of great interest since the bond MOs description has already produced very satisfactory 
estimations [10]. However, the methodology can be easily extended to more sophisticated 
Hamiltonians, including for instance one-site U and inter-site V repulsions, the resolutions of 
which are traditionally derived from expensive approaches. 
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